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, Abstract. We review the construction of homological evolutionary vector fields on infinite jet 

spaces and partial differential equations. We describe the applications of this concept in three 
(— I I tightly inter-related domains: the variational Poisson formalism (e.g., for equations of Korteweg- 

de Vries type), geometry of Liouville- type hyperbolic systems (including the 2D Toda chains), 
and Euler-Lagrange gauge theories (such as the Yang-Mills theories, gravity, or the Poisson 
sigma- models). Also, we formulate several open problems. 
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, Introduction 

' We present a systematic yet very basic review of the construction of homological evolutionary 

vector fields on the infinite jet spaces and on several natural classes of partial differential equa- 
tions such as the gauge models. For a long time, this geometric structure developed in parallel 
in mathematics and physics. Since mid-70's, it has been being used intensively in theoretical 
physics: specifically, within the BRST- or BV-technique for the quantization of gauge-invariant 
systems [71 l9l \25 \ \29 \ I59j. In mathematics, the concept first stemmed over smooth manifolds in 
the problem of the homological vector field realizations of Lie algebroids [6l] , which encompass 
Lie algebras and vector bundles and are best known in the framework of symplectic and Poisson 
' geometries. Although the notion of Lie algebroids in the geometry of manifolds had been known 

for a long time [491 , the proper generalization of this omnipresent structure to the geometry 
of jet bundles (c.f. [46) ) appeared in [37]. In principle, it could be contemporary to the discovery 
of the BRST-technique ([Dj vs [46]). In the meantime, a steady progress in the geometry of jet 
spaces [231 [US] created the platform for important applications of the future concept. Let us 
name only a few of them: 

• the variational Poisson dynamics for the KdV-type systems [I51|23], also [1U| : 

• modern revisions of the theory of Liouville- type hyperbolic equations [15U67j . in particular, 
the 2D Toda chains associated with the Dynkin diagrams [48] : 

• the cohomological approach to gauge fields [3l 0] and (quantum) Poisson theory [10] . 

In this paper, we recall the jet-bundle geometry that stands behind the structures = 0. It 
allows us to address, from a unified standpoint [SU [59], many relevant physical models ranging 



from the Korteweg-de Vries equation and its Poisson structures [50] to the Yang-Mills theories, 
gravity, and Poisson sigma-models [U [12] . 



three above-mentioned classes of geometries are regular sources of collections of such operators. 
We impose the Berends-Burgers-van Dam hypothesis [11]: namely, we let the images of the 
operators be subject to the collective commutation closure. This involutive setup gives rise to 
the bi-differential structure constants [19j and bi-differential analogs of Christoffel symbols [38j. 
For any number > 1 of the operators Ai, An, we reduce the setup at hand to the 
variational Lie algebroids |37j with one variational anchor. Having endowed the new bundle 
geometry with the Cartan connection, we reverse the parity of the fibres in it by brute force 
and then construct the odd, parity-reversing evolutionary vector field Q that encodes the entire 
initial geometry via the homological equality = 0. 

In section [3] we describe three natural examples of such geometries: the variational Poisson 
algebroids [271 13Z]) 2D Toda chains viewed as variational Lie algebroids [371 IMI, and gauge 
algebroids [2l[38]- 

First, in section 13.11 we recall the construction of necessary superbundles using the notion 
of the variational cotangent bundle [U] and derive the representation via = for the 
variational Poisson structures [23]. As usual, the renowned KdV equation [50] offers us the 
minimally possible nontrivial illustration. The odd vector fields Q are themselves Hamiltonian 
and we write the corresponding PF-charges explicitly. 

By exploiting the profound relation between the variational Poisson geometry and 2D Toda 
chains [15l |35l |39l , we obtain the evolutionary fields Q for these models of Liouville type [62] . 
We emphasize that the differentials Q and their cohomology carry much more information 
than the ordinary Chevalley-Eilenberg differentials for the semi-simple Lie algebras at hand. 
Being obtained from the gauge-invariant Yang-Mills equations under a symmetry reduction [48j , 
the 2D Toda chains represent the vast class of nonlinear Euler-Lagrange systems of Liouville 
type [581 ESI EZ]. This class is very interesting by itself; moreover, the construction of the 
homological vector fields Q in section [32] remains valid uniformly for all such equations. 

Finally, in section 13.31 we expose the true geometric nature of the BRST-differentials [9] or 
the 'longitudinal' components of the BV-differentials [^ for gauge systems. In agreement with 
the Second Noether Theorem, the differential operators Ai, . . ., An which we study here emerge 
from the Noether relations between the equations of motion; we explain why the same approach 
of Q2 = can successfully grasp more subtle geometries without any further modifications. 
The gauge parameters (that is, the arguments of the operators Ai), being the parity-reversed 
neighbours of the ghosts, together with the antifields (which will have appeared before in 
section 13. ip and the newly introduced antighosts constitute the classical BV-zoo. Similarly 
to the variational Poisson case, the evolutionary BV-differentials D, which incorporate the 
odd vector fields Q, are then determined by the Schouten bracket of the BV-charges with the 
inhabitants of the zoo. 

In the appendix, we report the results of a direct search for the linear differential operators 
with involutive images; their coordinate-free understanding as the variational anchors in 
variational Lie algebroids is yet to be achieved. 

We conclude that the presence of the evolutionary differentials Q is not specific to either the 
variational Poisson or gauge systems alone: indeed, it is an immanent property of a much wider 
class of models. The study of the arising cohomology theories not only sheds new light on each 
system or structure in particular but also reveals far-reaching links between them. 




1. Preliminaries: Notation and conventions 

This towel may not serve as a shelter during the hurricane. 

(Liability disclaimer) 

The geometry of jet fibre bundles is an established and well developed domain [3l HSJ [511 [531 EH 
I65j . Any deviations between the notations from different sources are easily tractable; we follow 
the notation of \37\ [38] and adopt Dirac's convention on bra- covectors ( | and -ket vectors | ), 
fixing the orientation 5p A 5q on the cotangent bundles. 

By default, we pronounce the standard incantation about the ground field M, the infinite 
smoothness of all mappings at hand, the locality along the oriented base of the bundles 
(whence emerge the jets of sections of vector bundles vr : — >. instead of the jet spaces 
joo^^n _^ M"^) for maps of manifolds), and the locality with respect to the jet order (whence 
the filtration in C°°( J°^(7r))). We assume that all total differential operators are local, i.e., 
polynomial in the total derivatives. By convention, the symmetries are infinitesimal and we 
identify the evolutionary vector fields | d^) with their generating sections | (p). For the sake of 
transparency, we consider the purely even, non-graded initial geometry of commutative vector 
spaces, manifolds, and bundles; the odd objects which we deal with appear later, through the 
explicit use of the parity reversion 11. 

Our reasonings are independent of the choice of local coordinates. We denote the independent 
variables by a; G whereas the notation q runs through all the dependent variables such as the 
gauge fields (if any in the model under study) at the points of space-time or say, the realizations 
of the "string" S"' in the space-time M^'^. The notation is multi- faceted: in each class of 
examples (see section [3D, we shall recall the common specific notation so that, on one hand, an 
object in this review may acquire several interpretations in different contexts, while on the other 
hand, we shall of course meet with the habitual ghosts, antighosts, and antifields. 

We work with the infinite jet bundles from the very beginning. Let us remember that many 
elementary constructions on usual manifolds -first of all, the Leibniz rule in the commutator [ , ] 
and the Schouten bracket | , ] or in the Poisson brackets under the multiplication of the 
arguments by smooth functions- either entirely disappear over the jet spaces or (seldom) 
survive after a careful consideration and proper amendments. Furthermore, we deal with the 
variational Poisson structures [23l[33] over the spaces J°°{Tr) of jets of sections for vector bundles 
vr: E^~^^ rather than suppose that the target space M"^ for J°°{E^ M™) alone is a 

Poisson manifold (c.f. [H [121 [22] ) . The variational Poisson theory we exploit does not amount 
to the canonical formalism which one applies to a Lagrangian system after taking its Legendre 
transformation (except for section [3^2] where we do use that approach as well), c.f. [5l [52]. 

For completeness, we explicitly endow the horizontal infinite jet bundles [33] over the 
spaces J°°(7r) with the Cartan connection. When introducing the modules of sections of the 
horizontal bundles such as the spaces of the gauge parameters, never we attempt to treat these 
sections as the "collections of functions" or use any other, equally weak interpretation of their 
structure because it does not grasp properly the geometry of the building blocks — instead, this 
shadows it with irrelevant fragmentations. Indeed, let us remember that, under a differential 
change of the coordinates x and q, the genuine reparametrization rules for those sections are 
imperatively prescribed by their geometric nature (see [371 [22] for detail and examples). We 
thus ought to be careful with the C°°( J°°(7r))-module structure for the spaces of such sections; 
this is specific in particular to the variational anchors in variational Lie algebroids (see section [2] 
and |37| ) which are neither the anchors nor Lie algebroids for usual manifolds, respectively, unless 
the base manifold E"' shows itself in space-time as no more than a point, or unless all gauge 
fields are uniformly constant. This also means that we deal ab initio with the field-dependent 
gauge parameters p[x,[q]y, moreover, their rigorous introduction (see section [3^3]) demonstrates 
that the gauge models must be particularly restrictive in order for these sections p to not depend 
on a part of the variables that encode the points of J°°(7r). 



In sections [2] and ISTTl we work on the empty jet spaces; later, in sections 13. 21 and l3.3t we do on- 
shell. We assume the off-sheh vahdity of the Berends-Burgers-van Dam hypothesis stating the 
collective commutation closure of the evolutionary fields "9^1^ whose generating sections Ai{-) 
belong to the images of linear differential operators Ai. This may require us to quotient out the 
trivial symmetries that vanish on-shell. 

We recall that the Euler-Lagrange systems £el = {Fi := 5S/5q^ ~ contain as many 
equations as there are unknowns in it. (In contrast with [21j, we do not assume that the density 
of the action functional 5 is a differential polynomial: indeed, see section 13.21 ) We point out 
also that usually the equations, either in evolutionary systems or in the Euler-Lagrange systems, 
are enumerated (more precisely, labelled) by the respective unknowns which explicitly occur 
in the left-hand sides. However, we emphasize that the admissible reparametrizations of the 
fields q and of the equations F = Q are entirely unrelated. This produces its due effect on 
the transcription of the infinitesimal symmetries q = Lp{x,[q\) in the former case and, in the 
latter, on the objects that lie in or are dual to the horizontal module Pq 3 Fi of the equations, 
e.g., on the antifields q^ (thus, more appropriately denoted by -F^), the generating sections -0 of 
conservation laws (for the Euler-Lagrange systems, emerging from the Noether symmetries ^pc 
by the First Noether Theorem), or the objects which obey the Second Noether Theorem such as 
the ghosts 7 and the antighosts 7^ . The 2D Toda chains in section 13.21 offer a perfect example 
of such discorrelation between ipc and which is brought in by force due to a purely esthe- 
tic tradition of writing hyperbolic systems with their symbols cast into diagonal shape. This 
requires a bit of attention (see [351 [32] ) : the Noether map — >■ ■0 is not always the identity. 

In this paper, we understand as the Noether maps not only the linear differential operators 
Pq 3 ip ^ if ^ sym.£ but, for i > 0, the operators Pi symf at each ith generation of the 
syzygies G Pi, i.e., the equations £ = {F = \ Fj € Pq} at i = 0, Noether's relations 

= at i = 1, etc.; the notation Pi corresponds to the modules of sections dual to Pi. In 
this sense, the Hamiltonian operators A: Pq ^ sym £ and the generators Ai: Pi sym £ of 
gauge symmetries are Noether operators. 

On the same footing, we note that the arguments Pi of the linear differential operators Ai 
which determine the evolutionary vector fields d^^ip ) possess an ambiguity. This is due to 
several independent reasons and occurs not only because of the on-shell equivalence: for a given 
equation £, we let p'^ ~ p'^ whenever p'-—p'lKiOon£. But indeed, the kernels ker Ai in reducible 
gauge theories carry the freedom as large as the presence of the free functional parameters in 
their description if the operator equations o Vj = admit nontrivial solutions Vj 7^ 0. 
Independently, the differential equations J2iLi Ai[pi[x, [q{x)fj^ = may have non-empty non- 
linear spaces of solutions q{x). This shows that the structure constants ([2]) (c.f. [I9] and [38j), 
which are bi-differential operators with respect to the arguments Pi, . . ., of the operators Ai, 
. . ., Aj^, are in fact equivalence classes of mappings. Because of this, considerable efforts are 
required to achieve the canonical normalization of the homological evolutionary vector fields on 
the jet-bundle extensions of such domains, see [37] and [2T] . 

The construction of gauge symmetries, which appear in a class of geometries very particular 
by themselves, is reviewed in section [3^31 But let us recall in advance that these symmetries have 
the form 5^^^ with arbitrary Pi = Pi(x, [q]). In this context, the class of hyperbolic partial 

differential equations q^y = f{x,y;q^,qi,qy) that admit arbitrary functions p{x) and p{y) in 
their symmetries is close to but different from the renowned class of Liouville-type hyperbolic 
quasilinear systems (see section 13.21 or [671 E2] and [33 [39] ) . 

Summarizing, we always highlight the geometric nature of the structures at hand and describe 
the links between them, as carefully as possible. The literature which is available on the 
rapidly developing topic of this review is so vast that, instead of the botany of new remarkable 
advancements, we now do the sociology of repeatedly practised oversights. 



2. Basic concept: Involutive distributions of operator-valued evolutionary vector 
fields 

Let vr and ^j, i = 1,. . . ,N, be vector bundles over the base manifold 3 x and denote by 
TToo : J'^i'^) and ^j^oo : J°°{Ci) the respective projections (the forgetful maps) to the 

base of the spaces of infinite jets of sections for these bundles. Consider the horizontal modules 
^i'^ooiCi)) of sections of the induced bundles 7r^(Ci) and denote their elements by Pi[x, [q]). 
The main example of such modules is x{tt) := r(7rSo(vr)) 3 if, which is the C°°(J°°(7r))-module 

of the sections of evolutionary derivations d!^^ ~ ^ '5^ ^ a^('/') af — " " " '^°°(^); we shall 

also use the modules A'^(7r) of the highest horizontal form^ C = L ■ 7r^(dvol (X"^) and the 

modules x(7r) of the variational covectors ip which are dua0 to the sections ip € x(7r) with 
respect to the coupling (, ) that takes values in A (vr). (Generally speaking, the construction 
of the bundles can be very involved even in the regular and natural classes of geometries, see 
section [32] and [39\ I37j.) Not accenting it at all times, we shall project to the nth horizontal 
cohomology spaces H (tt) that are determined by the lifting d onto J°°{tt) of the de Rham 
differential on U"'. By this we gain the right to integrate by parts; we denote by J w or by [uj] 
the equivalence classes of differential forms w. 

Suppose that Aii r(7r^(^j)) — )• x(7r) are linear total differential operators that take values 
in the space of evolutionary vector fields for each i = 1, . . . ,N. That is, the domains of the 
operators may be different but the target space is common to them all. The Hamiltonian 
operators A: x(7r) — )• >f(7r) on empty jet spaces J°°{7r) are an example of such structures [23]. 

Assume further that the images of the A*" operators are subject to the collective 
commutation closure, 

[im^i,imAj] C ^^^^im^^. (1) 

The sum of the images in the right-hand side is not direct because in effect they can overlap [50J . 
Assumption ([1]) gives rise to the bi-differential structure constants [19j via 



N 



[Ai{pi),Aj{pj)] =^Ak{c^j{pi,pj))- 



k=l 

each c^j is the equivalence class of bi-differential operators (possibly, with nonconstant, field- 
dependent coefficients) with both its arguments and its values taken modulo the kernels of 
the operators Ai, Aj and A^, respectively. (The kernels can be sufficiently large unless the 
extra nondegeneracy assumptions are made; for the sake of brevity, we shall not emphasize the 
presence of such kernels in all formulas but rather in few ones, c.f. |37].) 

Using the permutability of the evolutionary vector fields with the total derivatives and taking 
into account the Leibniz rule that always holds for the evolutionary derivations d^^^p ) ; we obtain 
the canonical decomposition of the bi-differential structure constants, 

4(P.,P,) = d^^l^^iP,) ■ - d^I^^p/p,) ■ 5f + r^(P.,P,). (2) 

The symbols F^^ absorb the bi-linear action of c^^- by total differential operators on the arguments; 
the notation is justified [38] because the terms F^^- do indeed transform as the Christoffel symbols 

^ The C°°(J°° (7r))-module structure for the spaces of Hamihonians H or Lagrangians £ is practicaUy useless 
because it does not pass -via any reasonable Leibniz rule- through the variational Poisson or Schouten brackets. 
^ The agreed choice of the coupling implicitly prescribes that we fix the metric and the volume form on the 
base E". Let us add that in this paper we never focus on the concrete choice of that metric and also we do not 
pay any particular attention to the upper or lower location of the indices. 



under the reparametrizations Pi ^ gPi invoked by a differential change of the jet coordinates 
(whence the g^s are hnear operators in total derivatives). We emphasize that generally speaking, 
the symbols T^j alone -that is, without the first two standard terms in ([2])- do not determine 
any Lie algebra structure: the true Jacobi identity is given by ([7]) where the standard terms are 
taken into account properly. It is likely that the restoration of the standard component in ([2]) 
for known bi-differential symbols T^j is the cause of a fierce struggle in the modern double field 
theory and in the theories which are based on Courant-like brackets. We notice further that the 
bi-differential Christoffel symbols introduced in ^ must not be confused with the connections 
A. = Ai dx* in principal fibre bundles in gauge theories. It is interesting to pursue further this 
emerging parallel between the geometry over the jet bundles and the known connection-based 
structures over usual manifolds. 

We remark that, in general, c^j ^ (mod ker A^) ii k ^ i, that is, the image of an operator 
alone may not be involutive even if the operator (in particular, its domain of definition) is well 
defined regardless of the entire collection Ai, . . ., A]\}. Also, we note that for only one operator 
{N = 1) the Lie algebra of evolutionary vector fields in its image is, generally, non-abelian, 
which must not be confused with the abelian gauge theories (e.g.. Maxwell's electrodynamics) 
corresponding to one-dimensional Lie groups, which are always commutative. 

We finally notice that the nature of the arguments Pj is firmly fixed by the reparametrization 
rules Pj ^ gPi- Therefore, the isolated components of the sections p^, which one may be tempted 
to treat as "functions" and by this fully neglect their geometry, may not be well-defined as true 
functions. We repeat that this primitivization discards much of the information about the setup 
and shadows the geometry which we have addressed so far. 

We now reduce, in a regular way, the setup of many operators and their individual domains 
to one operator on the large domain [37]. Namely, we take the Whitney sum over J°^(7r) of the 
induced bundles 7r^((^j), which in down-to-earth terms means that we compose the new "tall" 
sections 

N 

P = I P2 I • • • I Ptv) G r(0_^^^^^ 7r*^{^.)) =: Tn{Q, (3a) 

k=l 

and understand the sections p as the arguments of the new "wide" matrix operator 

A={Ai\A2\...\An). (3b) 
Obviously, property ([T|) implies that 

[im^,im^] C imA. (4) 

By definition, operator (|3bp is the variational anchor in the variational Lie algebroid [37]; its 
total space is the quotient of the sum over J°°(tt) of the initial domains for Ai by the kernel 
keryl = {p G rr2(^7r) | Z]j=i^i(Pi) = O}. We emphasize that variational Lie algebroids are 
not Lie algebroids and the variational anchors are not anchors as described by their traditional 
definitiorll [MllSl] for usual manifolds (here, for the target alone but not J°°(Z'" M^)). 

The operator A transfers the commutator [ , ] in the Lie algebra of evolutionary vector fields 
to the Lie algebra structure [, ]a on the quotient of its domain by the kernel [211 [37]. It is now 
entirely obvious that for a given collection Ai, . . ., Ajsf of the operators with a common domain 
(and not necessarily involutive images, but this does not matter) the induced bracket [ , ]a may 



The problem of integration of Lie algebroids to Lie groupoids, which is solved for manifolds (see [13]), is open 
in the variational case over the jet bundles. 



not determine a well defined Lie structure on a single copy of that common domain (c.f. [2j): 
Indeed, it is the concatenations p of the N elements from the same space for which the new 
bracket actually appears. 

In fact, it is not the horizontal modules r(7r^(^j)) of sections Pi{x, [q]) of the induced fibre 
bundles vr^(^j) which we really need for the representation of the geometry at hand in terms of 
the homological vector fields, but it is the product J°°{S,n) '■= J^iO «^°°(^) o'^^r Z"" of the 
infinite jet spaces 

for7rande = ©-Iiei- This means that we operate with the jet coordinates Pj.^ 

corresponding to the multi- indices r instead of the derivatives ^^Pj(a;, [g]), whence the linear 
differential operators Ai become linear vector-functions of the jet variables Pj.^. (In practice, 
this often means also that the number of the "fields" doubles at exactly this moment.) The 
following diagram endows the total space J°°{(,tt) with the Cartan connection Vc (c.f. |55j). 



J°°(0 xjjr. J-(7r) TZ^r: J°°(7r) 

(5) 



j°°(0 : . i^" . 



This justifies the application of the total derivatives to the variables p (including each p^ alone) 
and, from now on, allows us to consider on J°°(^7r) the evolutionary vector fields with their 
sections depending (non)linearly on the fields q and the variables Pj. 

The conversion of the sections Pj to jet variables creates the miracle: the kernel ker A of the 
operator A becomes a linear subspace in J°°{£,n) for each point of J°°(7r). This is a significant 
achievement because it allows us to operate with the quotient spaces, which we anticipated when 
the bi-differential structure constants c^j were introduced in ([2]). 

Next, we take the fibres of the vector bundles (hence, of ^i,oo) and reverse their parity, 
n: Pj <^ bi, the entire underlying jet space J°°{tt) remaining intact. This produces the 
horizontal infinite jet superbundle J°°(n^7r) := J°°(n^) x^-n J°°(7r) — )• J°°(7r) with odd fibres, 
the coordinates there being bi ^ for i = 1, . . . , N and |t| > 0. (In a special class of geometries, see 
section [3^ we shall recognize the variables bi as the ghosts, which are denoted usually by 7j.) 
The operators Ai tautologically extend to J°°(n^7r) and become fibrewise-linear functions in bi-r', 
the equivalence classes of the bi-differential Christoffel symbols V'j^j are represented by bi-linear 
functions on that superspace. 

Theorem. Whatever representatives in the equivalence classes of the bi-differential symbols F^j 
be taken, the odd-parity evolutionary vector field 

1 ^ 



is homological: 



N N 

k=l k=l 



The proof is obtained immediately by passing from many arguments of many operators to the 
tall sections of wide operators, see ([3]), and then literally repeating the first half of the proof of 
the main theorem in [B7\. However, we still prefer to prove the theorem straightforwardly; this 
will attest that the nature of the variables bi at different i's may be entirely uncorrelated. 



Proof. The halved anticommutator = of the odd vector field Q with itself is an 

evolutionary vector field. In this anticommutator, the coefficient of d/dq, which determines the 
coefficients of d/dq^ at all \a\ > 0, is equal to 

N N N 

j=i k,e=i i,j=i 



we now double the minuend -taking into account that fej and bj are odd, whence the minus sign 
occurs after their interchange- 



TV N 



i,j=l k=l 

by the definition of the symbols T^j. 

Second, let us consider the Jacobi identity for the Lie algebra of evolutionary vector fields 
with the generating sections belonging to the images of the operators Ai, . . ., A]\f viewed as the 
fibrewise-linear functions on the total space of the bundle J°°{(,tt) J°°(7r): 



O 
{imn) 



O k 
{imn) 

where we have relabelled the indices £ ^ A; in the second sum, 



= - E E^4"^Sp,)(^--(p-'^'")) +Er|(rL(Pm,p„),p^)}- (7) 

(imn) 

Thirdly, the velocity of each odd variable fo^ induced by the (for convenience, not halved) 
anticommutator [Q,Q] = 2Q^ is obtained as follows, 

1 ^ 



Ai(6i) (E ^"-^^ 2 E ^mni^mi bn),bi) — - ^ ^ r^„(6„ 

m,n m,n i^i m,n 

N 



We notice that the extra sum over the three cyclic permutations of each fixed set of the indices 
i, m, and n does not produce any change of the signs because the cyclic permutations (of the 
respective three odd b's) are even. Consequently, by taking the sum of all possible Jacobi 
identities ([7]), we conclude that the sought-for coefficient is equal to 



= E{---}=° {mod(p = {ipi,...,ipN)\^Ak{vk) = 0). 

i,m,n O k=l 
(imn) 

This completes the proof. □ 



We remark that the above construction of the homological vector fields Q in absence of the 
equations of motion but with an a priori given cohection Ai, . . ., Ai\i of generators ([6|) for the 
infinitesimal transformations of the model patterns upon the construction of nontrivial gauge 
theories with the zero action functional and the prescribed nontrivial gauge group, see j8l|60j. 

3. Examples: Three regular classes 
3.1. Variational Poisson algebroids 

The definition of Hamiltonian operators on jet spaces |23j, which are a class of linear differential 
operators in total derivatives, reads as follows \53\ 133] , Consider a total differential operator 

A: >c{tt) x(7r) that maps variational covectors to evolutionary fields. Let Tii, 7i2 G IT^^tt) 
be two Hamiltonian functionals, that is, the equivalence classes of two elements in the 
coupling module A (vr) of the highest horizontal forms. The operator A is Hamiltonian if 
the bracket { , }a defined by the formula 

is a Poisson structure, i.e., if it is bi- linear (which comes automatically), is skew-symmetric 
(whence A = —A^; the Z2-graded case is more subtle, c.f. [31,j), and satisfies the Jacobi identity. 
There are many formulations of convenient criteria for the verification of the latter (see \53\ |23] 
and also [411131) ). It is important that they reveal -or are based effectively on- the commutation 
closure ^ for the images of Hamiltonian operators. 

For a Hamiltonian operator thus defined on the empty jet space J°°{tt), a determined 
autonomous evolution equation is called Hamiltonian with respect to that given structure if 
the equation can be cast into the form 

f = {F = g,-^(^)=0|^€F"(vr)}. 

The left-hand side(s) of the equation £ belong to a suitable horizontal module Pq of 
form r(7r^(,^)) for some ^; usually, this module is identified with the module >f(vr) = r(7r^(7r)) 
of the velocities. This identification conveniently fills in several small gaps in the count of the 
base dimensions, the location of the Hamiltonian H for £ in the highest horizontal cohomology 
group but not in the preceding one, and the (in/ex) elusion of the time-derivatives of q in the 
module Pq of equations. However, we emphasize that the reparametrizations of the unknowns q 
and equations F = are not correlated in principle. Let us remember that the arguments p 
of the Hamiltonian operator for £ belong to the module Pq which is ( , )-dual to Pq. While the 
system £ remains in evolutionary transcription, the variational derivatives Sg/dq of conserved 
densities g for £ constitute a linear subspace of Pq, but this correspondence vanishes under more 
general reparametrizations of the equations. 

Reversing the parity of the sections p, we arrive at the following diagram: 

antifields — ► b G HPq 

n 

p € Pq " " F ^ Pq ^ >c{tt) < — usual identification. 

We remark that, under the labelling of the evolution equations £ = {q^ = r.-h.s.} by the fields q 
themselves, the antifields are then, on these grounds, usually denote^fl by q^ or q*. But in 

* Without and with regard to the underlying metric and the Hodge structure, respectively (c.f. 12 for an 
example); because we do not have any metric involved explicitly in the model, we entirely discard this distinction. 



view of the non-identical correspondence between the fields q and the left-hand sides F of the 
equations £ = {F = 0}, the proper notation for the antifields would be F^ (resp., F*). 
Not only the images of Hamiltonian operators are involutive so that 

[^(Pi),^(P2)] = A{[p^,p^]a), 

but the brackets [ , ]a are calculated effectively for each given A, see, e.g., [HI [31]. Therefore, we 
have all the data of section [2] ready at hand: taking the odd jet variables b = n(p), we obtain 
the odd evolutionary vector fields 

which encode the variational Poisson structures (primarily, the Jacobi identity, see above) in 
terms of the homological condition = 0. 

The homological vector fields Q for Hamiltonian operators were constructed in the recent 
paper [27] by taking the advantage of the known geometric interpretation for b as the parity- 
reversed variational covectors -but without references to variational Poisson algebroids. Namely, 
this was done by explicitly calculating the induced velocities of the sections p[x, [q]) for a priori 
specified velocities (p = A{p) of the unknowns q; we shall apply the same technique in section [HTBl 

Example. There are many examples when the Hamiltonian operators A are addressed in the 
context of the brackets [, ]a and bi-differential symbols T^A- instance, the notorious KdV 
equation (say, upon w{x,t)) yields ([67], see also the Appendix) 

(9) 

which is the minimally possible nontrivial illustration; we refer to [35] for the Boussinesq system, 
the Kaup-Boussinesq equation was considered in [36]; the Hamiltonian operators and the bi- 
differential Christoffel symbols for the Drinfeld-Sokolov equations associated with the root 
systems of rank two are contained in |40j . 

The calculation of variational Poisson cohomology groups determined by the evolutionary 
Poisson differentials Q is a rapidly developing topic, see [MKH]. We regret to notice the extent 
to which the progress in the cohomology theory for variational Poisson algebroids is retarded 
-by three decades at the least- in comparison with the BRST- and BV-cohomology technique 
for the gauge algebroids (see section [331) • 

The evolutionary vector field Q in ([8]) is itself Hamiltonian with respect to the even VF-charge 
Cl = —^{b, A{b)) and the canonical symplectic structure oj = ( g) ^^^^ stems from the volume 
form 5p A 5q in the fibres over points of in the variational cotangent bundle |47| to the jet 
space J°°(7r): 

O = d^"^- + d^''^ - 
^ 5U/5b^ -5U/5q 

We emphasize that this property of the field Q associated with a Hamiltonian operator A does 
not require the target M in the initial geometric setup J°°{X!"' M™) to be a Poisson manifold, 
c.f. [ll|12l[22]. 

The condition = \[Q,Q] = upon the evolutionary vector field Q is equivalent to 
the variational classical master equation |ri,ri] = upon the variational Poisson bi- vector 

ri = ^(6, A(6)) and the variational Schouten bracket |a;i,a;2] = (^'^u>iA 6 u>2^ (the odd- 
parity Poisson bracket or the antibracket, see [68j or |32j and references therein), which in the 



adopted notation and with Dirac's convention reads 




where the brackets [ ] mark the cohomology class after the volume form dvol Z"" is fixed in 
the coupling (, ). The composition of the antibrackelH relies on the arrangement of the jet 
(super)fibre variables in the coupling-dual pairs of opposite parity: here, we have the fields q 
and the antifields b = q"^ so that (up to the minus sign in fi = —ft which is due to Dirac's 
convention) 

Gauge models, which we address in section [331 provide a different geometric interpretation for 
the antifields and require the introduction of at least one generation of the ghost-antighost pairs. 

The definition of variational Poisson structures admits its further generalization to non- 
evolutionary systems £. While the bi-linearity of the bracket remains intact, its skew-symmetry 
imposes a nontrivial constraint upon the operator A: Pq ^ symS in terms of the linearization 
of the system £ at hand, whereas the Jacobi identity amounts to the same condition = 
upon the odd vector fields dS]). A significant progress in this direction has been achieved only 
very recently, see [Ml US] • 

Open problem 3.1. Do there exist Hamiltonian operators Ai, . . ., A^ which are compatible 
in the sense of ([1]), and what then could be the geometry standing behind the bi-differential 
Christoffel symbols (c.f. [38])? 

Finally, we recall that the root systems of semi-simple complex Lie algebras are a regular 
source of variational Poisson structures for KdV-type models |I5]; this is achieved by using the 
2D Toda chains. 

3.2. 2D Toda chains 

The hyperbolic two-dimensional non-periodic Toda chains j48| [62] associated with the root 
systems of semi-simple complex Lie algebras g of rank r are the exponential-nonlinear Euler- 
Lagrange systems 

dj, = exp(if;g^), l<i<r, 

where K is the Cartan matrix of q. The action functionals for the 2D Toda chains are expressed 
in terms of the roots and are non-polynomial, thus exceeding the assumptions of [21]. Also, it 
is very instructive to inspect the non-identical correlation between the Noether symmetries and 
the generating functions of conservation laws for these Lagrangian equations, see |35[ [39] for 
detail. 

The 2D Toda chains are the best representatives of the vast class of Euler-Lagrange systems 
of Liouville type (see [671 [63] and references therein, also [351 [39]). By definition, the systems 
of this type possess as many Liouville's integrals w{x,y\ [q]) such that -^{w) ~ on-shell and 

similarly, w[x,y; [q]) satisfying ^{w) ~ 0, as there are unknowns q. For example, consider the 
Liouville equation £" = = q^y — exp(2g) = O} and let w = — qxx so that g^(^) = ~^(-^)- 

^ Let us recall that the Schouten bracket of variational one-vectors {b,ipi) calculates the commutator [(^1,(^2] = 
dip'i {'P2) — 01^2 (ifii) of the two evolutionary vector fields djpl' : Namely, we have that {{b, cpi) , {b, ip2)} — {b,[ipi, 1/^2] )■ 
This formula illustrates our earlier warning about the absence of any reasonable Leibniz rule for either [ , ] or | , ]. 



We have shown in [35j that the integrals for this class of Euler-Lagrange equations in fact 
depend differentially on the momenta m which emerge from the kinetic part of the action; this 
allows us to cast 8 into the evolutionary form qy = (5H/(5m, = —5Yi/5q. The differential orders 
of the integrals w = w\m\ the concide with the exponents of the Lie algebra q, see [Ml EH]; there 
are constructive procedures to obtain the integrals w and w (see [58] and [61] illustrated in [30] ) . 

In [35t [39] we proved that the adjoint linearizations of the integrals with respect to the 
momenta, 

U={t))^ (10) 

are linear differential operators with involutive images, see (H]). These operators produce 
symmetries of the Liouville-type system at hand (in particular, its Noether symmetries (pc = 
5'H\_w[m\\ / 5m. = \Z\{5'H[w]/5w)). By construction, each column of the operator □ stems from 
the respective integral for £. However, the image of a particular column may not itself be closed 
under commutation: we encounter an example for the root system A2 (here, for the second, 
higher-order integral and, respectively, the second column in ()10p . see [351 140j ). 

The geometric nature of the domains of the operators □ is the most nontrivial in comparison 
with the other two sample geometries which we discuss in sections 13.11 and 13.31 of this review. 
In brief, the arguments p{x, ['u;[m]]) in each such domain belong to the variational cotangent 
bundle to an equation other than £, while the Miura substitutions w = w\m\ relate the systems; 
we refer to [35] and [39] for detail. This reasoning implies that the components of the sections p 
do not exist individually as well-defined "functions" and therefore the operators □ may not be 
split in separate columns. (Let us remember that the operators whose domains appear from the 
variational cotangent bundles [37]; like the Hamiltonian operators, were referred to in |39[ I37j 
as the operators of second kind.) 

The term "Miura" is indeed justified in this context [15]: the brackets [, ]n on the domains 
of □ are calculated by using the correlation [33] of the entire geometry of the 2D Toda chains 
with the previous case of variational Poisson algebroids (section l3.ip . The correspondence yields 
the symbols explicitly in terms of the integrals w and the roots of g; this calculates the 
commutation relations in the symmetry algebras for the 2D Toda chains. 

As usual, we reverse the parity of the arguments p for the operators □ by taking 11 : p ^ 6 
and obtain the homological evolutionary vector field 

Example. The operator □ = g^, + for the Liouville equation qxy = exp(2g) specifies the 
differential Q = + 5^^^; the equality of the even velocity hbx of the odd variable h to 

the respective velocity which we obtained in ([9]) for the Korteweg-de Vries equation is no 
coincidence [39] . 

Open problem 3.2. Is it true that for each operator □ of second kind and for the associated 
differential Q there is the Hamiltonian operator A with the same domain such that = ? 

Finally, we recall that the 2D Toda chains related to the root systems are obtained by 
imposing the cylindric symmetry reduction in the Yang-Mills equations, the distinguished gauge 
models [48] . In turn, the latter do carry their own homological vector fields; we discuss them in 
the next section. However, the Yang-Mills equations are a class of models for which (unlike, e.g., 
gravity) the bi-differential Christoffel symbols standing in the fields Q are structurally simple. 
This shows the "preservation of complexity" in its genuine sense of logic under the reduction 
(large, simple) 1— ?> (small, complex) of the Yang-Mills systems to the 2D Toda chains. It is 
interesting to track the details of the correspondence between the homological vector fields Q, 
the compositions of the Schouten bracket, and the respective syzygies and charges in the two 



models and by this, to approach the origin of Liouville's integrals w along a new direction. This 
will give us a clue to their analogs in the afhne case of g^*^^. 



3.3. Gauge systems 

The Euler-Lagrange systems £ = {F = 6S/6q = \ S £ H (it)} contain as many equations 
as there are unknowns q. By construction, such equations Fi = are conveniently labelled 
by the respective fields q^; let us remember that the horizontal module Pq of the sections F is 
then Pq ~ x(7r). Because of this, the generating sections ^ (z Pq of conservation laws for £ 
acquire the nature of Noether symmetries (pc G sym£, which is indeed the case by virtue of 
the First Noether Theorem (up to, possibly, the non-identical Noether maps Pq >^('?i"), as we 
recalled in section I3.2p . 

Although the systems £ are determined, there may appear the differential constraints (also 
called syzygies, Noether identities, or Bianchi identities) between the equations of motion, 

$[F]=0 Vq = q(x), G Pi. 

The relations <I>i+i[<&j] = between the relations, valid identically for all ^i-i (here = F 
Pq) give rise to possibly several but finitely many generations of the horizontal modules Pi B $i 
for i > 0. In the sequel, we shall assume that the given action functional S determines the 
system £ = {F = 0} of equations of motion with one generation of the constraints ^[F] = and 
that there are no further relations between the already known ones. In the following diagram 
we summarize the notatiorH and interpret the objects at hand in terms of the BV-theory [7j: 



ghosts 



gauge parameters 



7, = 6i e nPi 



ei = Pi e Pi 



jI e Pi - 

id 

^'[F] £ Pi 
G nPo ^ 

. 

n 



Noether symmetries — >- ipc = £ x{'k) ~ Pq 



FGPn 



antighosts 

Noether identities 
antifields 

usual identification. 



Each relation <I>*[i^] = between the equations Fj = 0, 1 < j < m, yields the linear total 
differential operator Ai: Pi — )• symiS C x(7r) that generates symmetries of the system £. 
These symmetries are parametrized by arbitrary sections Pi{x, [qfj € Pi and are called gauge 
symmetries. 

Namely, suppose that the identity <I>[i^[q]] = holds irrespective of a section q = q{x). 
Consequently, this identity is indifferent to arbitrary infinitesimal shifts, which are given by the 
TToo-vertical evolutionary vector fields d\^'^ on J°^(7r). The chain rule implies that 



di'\F) 



($) = G Pi. 



Passing to the linearizations (the definition of which in the non-graded case is ^^■'((/^[a]) = 
d!^\tlj) whenever this right-hand side is well-defined for the given a and ijj), we conclude that 

^ Again, we do not fix any metric on S" and therefore not grasp the subtle difference between the transcripts 7^ 
and 7* for the antighosts (see, e.g., [12)1. We repeat that we do not pay any particular attention to the upper or 
lower location of indices. 



where the composition is a hnear mapping from x{tt) to Pi. Let us now couple this zero velocity 
along Pi with any element p € Pi from the dual module and then integrate by parts, staying 
in the equivalence class of the zero in the cohomology. We obtain 

(4'V°(4''V(p) = o. 

But let us recall that the linearization 

4) : x(7r) ^ Po ~ ^(^ 



and the adjoint linearization 

(4 V : ~ x(7r) ^ ^(^ ~ Po 

coincide for the Euler-Lagrange systems S = {-F = O} due to the Helmholz criterion 
(see [Ml [65]) 

= (^J))^ ^ 3 S £lt{Tr) \ F = 6S/6q. 

This implies that (^^^)\p) is, on-shell, a symmetry of the system £ for any section p € Pi. 
Moreover, for each Noether identity S Pi we have constructed explicitly the operator 

that yields the symmetries (p = Ai{-) of the model. 

We note that the constraints ^ need not be linear. Furthermore, in the cases of a very specific 
geometry of £ there may appear the symmetry-producing operators which do not issue from any 
differential relations between the equations of motion. The approach which we have formulated 
by so far treats all such structures in a uniform way. 

However, under the additional assumption that the symmetries obtained in the images of 
the operators are Noether, 9^(p)(>S') = [0] G 'H^{tt), the existence of the respective sections $ 
is justified easily and besides, these constraints appear to be linear. Actually, suppose that 
d'-f^p^iS) = {6S/6q,A{p)) = (p,At(F)> = [0] for aU sections p £ Pi. This implies that the 

Noether identity ^[F] = amounts to the linear differential relation A^{F) = between the 
equations of motion so that $ coincides with its own linearization. 

Example. The Maxwell equations = |P* = djj-^^ = O} upon the skew-symmetric field 
strength tensor Tij = diAj — djAi satisfy the obvious relation ^> = 5jP* = (that is, 
^ = {didjJ-^^ = 0}), which manifests the invariance of the system (including its action 
functional) under (finite or infinitesimal) gauge transformations A. f-^ A. — d^R^pix , [A])) of 
the electromagnetic field A = Aidx'^, here ddR is the de Rham differential on the Minkowski 
space-time. The conventional plus sign in the transcript of the Bianchi identities yields the 
minus sign in front of the adjoint operator d^j^ = — d^R in the gauge symmetries. Because the 
coefficients of the de Rham differential, which spreads along the diagonal of the (4 x 4)-matrix 
operator A, are constant, the bi-differential Christoffel symbols vanish in this model. 

It is remarkable that the Yang-Mills theories with non-abelian gauge groups still do constitute 
a regular class of examples with the field-independent (and at most constant, whenever nonzero) 
coefficients in the bi-differential Christoffel symbols r^^ (-, •). On the other hand, gravity produces 
the drastically more involved structure of these symbols with the explicit dependence on the 
unknown fields in their coefficients [20]. This is immanent also to the Liouville-type systems 
(see section 13. 2p . 



Let the operators Ai, Aj\[: Pi — sym£^ C >c{-k) be the entire collection of the gauge 
symmetry generators for the system £. For the sake of clarity, we always assume the off -shell 
validity of the Berends-Burgers-van Dam hypothesis pT] about the collective commutation 
closure for the images of these operators, see ([1]). Then the odd evolutionary vector field ([6]) 
is the Becchi-Rouet-Stora-Tyutin differential [9j; in the model at hand, all parity-reversed 
arguments belong to one module IIPi and are known as the ghosts, usually denoted by 7j. 
The variational Lie algebroid encoded by the evolutionary differential Q was named the gauge 
algebroid in 

As in the variational Poisson case (see section 13. ip . we arrange the variables in pairs of 
opposite parity: even fields q < — > odd antifields , and odd ghosts 7j < — > even antighosts 7^. 

The evolutionary BRST-field Q lifts to the odd derivation Q on the horizontal infinite jet 
superspace that now contains the antifield coordinates qX and similarly, the antighosts 7|.^, 
here |r| > 0. The induced velocities of the newly incorporated fibre variables are calculated in 
agreement with their geometric nature. 

First, we recall that the off-shell determining condition d!^\F) = V(-F) for ip to be an 
infinitesimal symmetry of the equation F = states the existence of the linear total differential 
operator V such that the right-hand side of the condition vanishes on-shell (i.e., by virtue of the 
equation and its differential consequences). Let us remember also that the operator V = 
depends linearly on ip because the entire equality does so. Specifically for the gauge geometry at 
hand, this means that d^^^^^ (F) = V^(^) (F) or, after the parity reversion U: F ^ F^^ = q^ of 

the linear entry F in both sides of the equality, we obtain the odd-parity component d^*^ ^ ... 

VA(^)(qT) 

of the lifting Q for Q. 

Second, we see that -due to the chain rule- the odd velocity of a Noether relation ^[F] G Pi, 
or of the even antighost 7^ associated with it, is equal to 

Corollary. The proper antifield -antighost lifting Q of the BRST -differential Q reads 

k ' > k /^v^AA-iAy^kiysq) + ^k\^ Adni)y&q>) j 

i ^ ' 

The full Batalin-Vilkovisky differential F) = dx-T + Q + ■ ■ ■ merges the parity-reversing (and 
besides, vanishing on-shell) Koszul-Tate differential. 



k 



with the lifting Q of Q; however, D may contain further terms of the correction that ensures 
the off -shell equality = 0. 

Let us now look at the problem of construction of the BV-differential from another 
perspective. The pairwise arrangement q, 7 -f-)- 7^ for the ( , )-dual variables of op|)osite 
parity prescribes the composition of the variational Schouten bracket |cji,a^2l = {6 ui A 6 0^2) 



( [68] and, e.g., [281 E]), ^ow extending to all generations of the ghost-antighost pairs: |a;i,a;2l 



6q^ 6q 6q ' J \ S-f ' 6^'^ ' 6j 



+ 



(^7 (57t ^7 



In this notationj^ the BV-differential D is determined (|68l [66], see also section I3.ip by the 
variational Schouten bracket of the BV-action S and the respective super-jet variables, 

where the extended action functional S = —S has even parity and S equals [181 [301 US], with 
the conventional Einstein's summation over repeated indices and Dirac's ordering of (co)vectors, 

5 = 5 + (q^ Afc(7fc)) - ^ + (correction terms) . 

The possible necessity to introduce the correction terms of higher polynomial orders in 7, 7^, 
or q^ is legitimate, e.g., when the coefficients of the bi-differential symbols T^j depend explicitly 
on the fields q and hence produce the redundant terms in the velocity of g^, which thus must be 
cancelled out properly. The renowned paper [12j (c.f. [22j) contains a perfect example of a gauge 
model with one generator A of its gauge symmetries and an explicit calculation of the (zero-order) 
bi-differential Christoffel symbols T^^, BRST- and BV-differentials, and the master-action S, 
see [12. p. 595] where the parity reversion 11 is the 'promotion' and the BRST-field Q is denoted 
by Sq. 

The cohomological condition [D,D] = is equivalent to the classical master equation 
fS, S] = 0. The BV-cohomology with respect to the differential D (or [S, •]) provides the 
resolvent for the algebra C°^{£) of observables in the model and opens a way for the quantization 
of gauge systems. 

Open problem 3.3. Describe the variational Poisson structures that are pertinent to gauge- 
invariant systems. 

The most recent techniques of [26^ [Ml HI] allow -or soon will allow- the regular search 
for variational Poisson structures on non-evolutionary gauge-invariant models by using the 
respresentation of the structures via = 0. It must be noted however, that the gauge algebroids 
are endowed a priori with the evolutionary differentials which suit well for the quantization of 
the systems. Therefore, the variational Poisson structures -whenever found for a given gauge 
model- will either provide only a few of its symmetries or, if there is a functional freedom in 
the domains of the Hamiltonian operators, will definitely not provide the entire Lie algebra of 
infinitesimal gauge transformations because of the apparent difference between the domains of 
definition for the Hamiltonian operators and ordinary gauge generators (compare the diagrams 
in this section and in section 13. ip . This confirms that the bi-Hamiltonian complete integrability 
paradigm ([50], see [S]), which relies on the construction of the Poisson pencils for hierarchies 
of evolutionary systems, becomes insufficient in the geometry of gauge fields. 

^ Likewise, the odd Laplacian is 

^ V ^ V 



■ o ■ 



5q^ 5q 5y Sf'' 



Conclusion 

The realization of nonlinear systems, or parts of the information which they carry, in terms of 
homological evolutionary vector fields Q is not specific to gauge-invariant models only. Indeed, 
other natural classes of partial differential equations can be addressed with convenience from the 
same viewpoint. Moreover, even for the already known gauge systems the available approach is 
ready to capture more of their geometry without any extra modifications. 

For a long time, the physical and mathematical incarnations of the notion of Lie algebroids 
over infinite jet spaces developed in parallel, the progress on the math side being seriously 
retarded with respect to the demands from physics. Seemingly, the exchange between the two 
theories, foresighted in pil], channeled through the construction of the variational Schouten 
bracket. Nowadays, the systematic and well-motivated calculation of the variational Poisson 
cohomology only begins — while the BRST- and BV-cohomology theories are an established 
domain of research and the subject of a vastest literature. In this review we traced the links 
between the variational Poisson algebroids, gauge algebroids, and -playing the role of mediators- 
2D Toda-like systems also viewed as variational Lie algebroids. 
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Appendix: The hunting of the variational Lie algebroids 

We now list several operators which we obtained in a fixed system of local coordinates and 
whose images are then closed under commutation. We denote by u{x) the field and, for brevity, 
we use the notation {{p, (7}}^ instead of T^j^{p,q). Let us fix the weights \u\ = 2, |d/da;| = 1 
that originate from the scaling invariance of the KdV equation ut = —\uxxx + 'iuux', we have 
that |d/dt| = 3. Using the method of undetermined coefficients, we performed the search 
for scalar operators that satisfy (jH) and which are homogeneous with respect to the weights 
not greater than 7. We obtained two compatible Hamiltonian operators A^'^^ = d/dx and 
^Kdv ^ _i(d/dx)3 + 2ud/dx + Ux, the odd powers (d/da:)^"""^^ of d/dx, and the Hamiltonian 
operator 

(as) + 3"^^ (^) + ^ + UUxxx- 

Also, there are four non-skew-adjoint operators with involutive images, 
- ul, Af^ = 2ul - uuxx - ^uux + -u^ (^)^ ' 

"^Ux ■ (pQx - PxQ), {{P, q}} a(6) = -2Ux ■ (pQx - PxQ) + U ■ {pQxx " PxxQ); 

^5 

^ ~ ^ - {-^) - Au^ {-^) , {{p, q}}^m = ■ {pQx - PxQ); 

-2uxUxx -ul-^, {{p, q}}^(7) = 8Uxx ■ {PQx - PxQ) + '^Ux ■ (pQxx - PxxQ)- 



A 



(6) 



{{p,<l}}. 



A. 



(6) 

(7) 



A, 



(7) 



Finally, we have found the operators that contain arbitrary functions: /(ti)(d/dx)" and f{u)v? 
with vanishing Christoffel symbols {{ , and also 

^3 = f{u)u^, {{p,9}}a3 = f{u){pxq-pqx); 

A4 = f{u) {u-^- Ux) , {{p, q}}Ai = f{u) {pQx - Pxq) ■ 

The transformation rules for the sections that constitute the domains of the operators A^, A4, 
A^\ . . ., Ag^ are yet to be found. Likewise, the interpretation as variational anchors and the 
understanding of the physics of the models remain an open problem for the infinite class of linear 
total differential operators with involutive images which were found in [56\ 167] , 
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